Bound states of a quartic and sextic inverse-power-law potential for all angular momenta

We use the tridiagonal representation approach to solve radial Schrödinger equation for an inverse-power-law potential of a combined quartic and sextic degrees all angular momenta. The amplitude singularity is larger than that but signs are negative positive, respectively. It turns out system has finite number bound states, which determined by ratio two amplitudes. solution written as series square integrable functions in terms Bessel polynomial.